Spin liquid behaviour in Jeff = 1/2 triangular lattice Ba3lrTi209 
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Abstract 

Ba3lrTi209 crystallizes in a hexagonal structure consisting of a layered triangular arrang ement of Ir^+(Je// = 1/2). Mag- 
netic susceptibility and heat capacity data show no magnetic ordering down to 0.35 K inspite of a strong magnetic coupling 
as evidenced by a large Curie- Weiss temperature 9cw ~ — 130K. The magnetic heat capacity follows a power law at low tem- 
perature. Our measurements suggest that Ba3lrTi209 is a 5d, Ir-based (Je// = 1/2), quantum spin liquid on a 2D triangular 
lattice. 



PACS numbers: 75.40.Cx, 75.45.+j, 75.47.Lx 

Introduction: Since Anderson proposed the resonat- 
ing valence bond model [l| , researchers have been search- 
ing for experimental realizations of quantum spin liquids 
(QSL) in geometrically frustrated magnets. In such 
materials, incompatibility of local interactions, called 

] frustration, leads to a strong enhancement of quantum 
spin fluctuations and effectively suppresses the long range 

] magnetic ordering. As a result, the material remains 
paramagnetic down to very low temperature compared 
to the Curie- Weiss (CW) temperature 0cw- The frustra- 
tion in these materials often arises from some special ge- 
ometries like triangular, kagome, pyrochlore, garnet etc. 

: i- 

\ The spin liquid candidates found so far have been 
mostly 3d transition metal based materials. A few 
examples are two-dimensional (2D) Kagome systems 

: SrCrgpGai. 9pOi9 {S = 3/2) [|, and ZnCu3(OH)6Cl2 
(5 = 1/2) 0, 5* = 1 2D triangular lattice antiferromag- 
net NiGa2S4 @, organic materials like 5* = 1/2 trian- 

■ gular lattice k-(ET)2Cu2(CN)3 etc. There are very 

\ few examples of spin liquid systems with Ad or bd spins. 
Na4lr308 a S' = 1/2 spin liquid in a three-dimensional 

] (3D) hyperkagome network, is probably the most notable 
member of the bd spin liquid family. 

Recently, Ba3CuSb209 {S = 1/2) with hexagonal 
space group P63mc was suggested to be in the QSL 
ground state High pressure hexagonal (PGsmc, 6H- 
B) and cubic (Fm-3m, 3C) phases of Ba3NiSb2 09 have 
\ also been suggested to be in the 2D and 3D QSL ground 
state, respectively (lo| . 

We have been searching for QSL candidates among 
hexagonal oxides with Ad/hd elements instead of 3cZ el- 
ements. The bd materials are very different from 3c? 
materials and thus interesting because of a weak onsite 
Coulomb energy but a strong spin-orbit coupling. For ex- 
ample, Sr2lr04 [llj and Ba2lr04 [l3| are reported to be 
spin-orbit driven Mott insulators. The magnetic prop- 
erties of these systems have been described based on a 
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Jeff = 1/2 state for the Ir^"*" ion. Among the various Ir- 
based compounds, Ba3lrTi209 is rather interesting since 
it has a chemical formula similar to the Cu and Ni-based 
compounds (discussed in the previous paragraph) and 
it crystallizes in a hexagonal structure [l^. However, 
detailed structural parameters have not been reported. 
Bryne et al. [l3| reported magnetic susceptibility of 
Ba3lrTi209 in the temperature range 77 — 600 K. High 
antiferromagnetic Weiss temperature (| 0cw |> 400 K) 
was obtained by them, suggesting that the magnetic Ir^+ 
ions are strongly coupled with each other. An obvious 
question arises, do they order at lower temperatures? If 
not then, is it a spin liquid system and a bd analog of 
Ba3CuSb209? 

Here we report preparation, structural analysis, mag- 
netic susceptibility and specific heat measurements on 
Ba3lrTi209. It crystallizes in space group PGsmc like 
Ba3CuSb209 and 6H-B phase of Ba3NiSb209. A large 
negative ^cw is obtained from CW fitting of susceptibil- 
ity data but no magnetic ordering is found from suscep- 
tibility and heat capacity measurements down to 0.35 K. 
Magnetic heat capacity follows a power law at low tem- 
perature. This indicates that the system is highly frus- 
trated and an example of a 5d QSL. We suggest that this 
is the first candidate of a bd based quantum spin liquid 
on a 2D triangular lattice with J^ff = 1/2. 

Experimental Details: Polycrystalline sample of 
Ba3lrTi209 was prepared by conventional solid state re- 
action method using high purity (99.9%) starting mate- 
rials. 

Powder x-ray diffraction (XRD) measurements were 
performed at room temperature with Cu radiation 
(A = 1.54182 A) in a PANalytical X'Pert PRO diffrac- 
tometer. Magnetization measurements were performed 
in a Quantum Design SQUID Vibrating Sample Magne- 
tometer (SVSM). Heat capacity measurements were car- 
ried out in the temperature range 0.35 — 295 K and field 
range — 9 T in a Quantum Design Physical Properties 
Measurement System (PPMS). High temperature (upto 
800 K) susceptibility was measured using a PPMS VSM. 

Results and Discussions: XRD measurement was 
done to check the phase purity of the sample and to de- 
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Figure 1: (i) Structure of Ba3lrTi209 without any site dis- 
order between Ti*^ and Ir^+. The triangular arrangement 
of Ir*"*" spins in the ab plane is shown (ii) IrTiOg dimers are 
shown (iii) One possible arrang ement of Ti*+ and Ir*+ ions m 
the ab plane is shown when about 1/3 of Ir"*"*" ions from Ir(l) 
site are exchanged with Ti''"'" ions of the Ti(2) site. 



termine crystal parameters, as t he p arameters were not 
mentioned in the earlier report |13| . The Ru-analog of 
Ba3lrTi20g i.e.. Ba3RuTi209 has been mentioned in lit- 
erature and it crystaUizes in the hexagonal P63 mc space 
group On the other hand, with a different Ir and Ti 
ratio, Ba3TiIr209 has been suggested to crystallize in the 
space group P63 /mmc [l5| . In both these space groups 
metal-metal structural dimers (2b sites or 4f site) are sep- 
arated by the 2a site metal plane. In PG^mc, the metal 
ions within the dimer are ordered while in P63 /mmc 
space group the metal ions within the dimers are not or- 
dered. We tried to refine our XRD data using these space 
groups and found P63 mc gives better refinement with a 
large site sharing by the Ti^+ and Ir'*+ ions (see sup- 
plemental material The lattice constants obtained 
from refinement are a = h = 5.7214A and c = 14.072lA. 

In the ideal case (i.e., without any site disorder), the 
Ti(3) site is occupied by Ti^+ ions and the Ti(2) and Ir(l) 
sites are occupied by distinct metal ions Ti^+ and Ir'*+ 
ions, respectively. This is indeed (nearly) the situation in 
Ba3CuSb209 where the Cu site is occupied only by Cu^+ 
(leaving aside a 5% site disorder), and Sb^+ ions are lo- 
cated at Sb(l) and Sb(2) sites. However in our case, we 
found a (37±10)% site sharing of Ir^+ ions with Ti^+ ions 
between Ir(l) and Ti(2) sites and (7 ± 4)% site sharing 
with Ti'*"'' ions in Ti(3) site. This is in fact not unex- 
pected, as their ionic radii are very similar. Sakamoto et 
al. also found 21% site sharing between Ti'*+ and Ir^+ 
in Ba3TiIr209 [l5| and a similar site disordered situa- 
tion is reported in the case of Ba3RuTi209 by Radtke 
et al. [131 • They studied probability of different Ru'*+ 
and Ti^"*" combinations based on high resolution electron 
energy loss spectroscopy and first principles band struc- 
ture calculations and concluded that site sharing of ions 
in 2b sites (i.e., Ir(l) and Ti(2) sites in our case) is more 
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Figure 2: Left axis; Magnetic susceptibihties of as prepared 
(AP) and quenched (QN) Ba3lrTi2 09 samples are shown. 
Solid lines denote fitting with CW law in high temperature 
range > 150 K. Right axis: Inverse susceptibilities (after sub- 
tracting xo) as function of temperature for both AP and QN. 
The dashed line is a linear extrapolation of high temperature 
data of AP sample. 



probable while site sharing with 2a sites (i.e., Ir(l) and 
Ti(3) sites in our case) is less probable. This was sug- 
gested because structural dimers of like ions (Ti-Ti) are 
energetically unfavourable due to a strong Ti-Ti repul- 
sion. The same reason is probably valid in our case and 
results in a small 7% site sharing between Ir^"'" ions at 
Ir(l) site and Ti'^+ ions at Ti(3) site. 

In case of perfect ordering among Ti'*+ and these 
two ions form face-sharing IrTiOg bioctahedra (shown in 
Fig. [TJii)) and Ir'*+ spins form an edge-shared triangular 
lattice in the ah plane, as shown in Fig. [TJ As a conse- 
quence of site disorder, the edge-shared triangular planes 
will be depleted. Further, Ir occupying the Ti(2) sites 
might also form a depleted triangular plane. A possi- 
ble arrangement is shown in Fig. [TJiii). The blue atoms 
represent Ti and the red atoms are Ir. 

Zero field cooled (ZFC) and field cooled (FC) mag- 
netic susceptibility was measured with different fields in 
the temperature range 2 — 400 K. No magnetic ordering 
is found down to 2 K but with 100 Oe field ZFC-FC split- 
ting is seen below 80 K (shown in Fig. [8] in supplemental 
material [lB|)- However, the splitting is very small (only 
11% of total magnetization at 2K) and supressed when 
measured even with 500 Oe. This suggests that a small 
fraction of the spins take part in a glassy state while 
the majority of the spins do not. In Na4lr308 also a 
small ZFC-FC (< 10% of total magnetization) splitting 
was observed below 6 K which the authors concluded as 
coming from a small fraction of the spins Fig. [2] 
shows the temperature (T) dependence of dc magnetic 
susceptibility of the as-prepared sample (light-blue open 
squares). Data obtained with field 5kOe using a SVSM 
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(2 — 300 K) and using a VSM with a high-temperature at- 
tachment with field 50 kOe (300—800 K) have been shown 
together. Susceptibihty data can be fitted weU with the 
CW formula in the high temperature (150 — 800 K) region 
(shown in Fig. (2]), which yields temperature independent 
susceptibility xo — 0.61 x 10"'^ cm^/mol, Curie constant 
C = 0.149 cm^ K/mol and Ocw = -133 K. In many Ir 
based oxides vn is found to be large and varies within a 
wide range (isl . [T8| . The C value obtained from fitting 
is much less than that expected for S = 1/2 magnetic 
moments (0.375 cm"^ K/mol) value. The large 9cw value 
suggests that there still are significant correlations in the 
triangular planes despite the depletion. The suppression 
of magnetic moments could be an effect of the extended 
nature of the 5d orbitals and the strong spin-orbit cou- 
pling expected for 5d transition metal oxides. Indeed, in 
the magnetically ordered iridates such as Sr2lr04 the low 
temperature saturation moments have been found to be 
less than one tenth of a /iB and effective moment in the 
paramagnetic region is found to be ^ 0.4/iB [l9| . 

With the aim of investigating as to how the prepara- 
tion procedure might affect the site ordering and hence 
the magnetic properties, we quenched the as-prepared 
sample in liquid nitrogen from 1000*^C. Comparing the 
normalised x-ray diffraction pattern of both as-prepared 
(AP) and quenched (QN) samples, we found that width 
of all the peaks and peak height of many peaks are de- 
creased in the QN sample. This indicates that the crystal 
symmetry is unchanged but ionic disorder (and possi- 
bly distortions) are less in the QN sample compared to 
the AP sample. Refinement of XRD pattern is consis- 
tent with a ^ 35% site disorder between Ir^+ and Ti'^^ 
cations at the 2b site but without any site sharing with 
the 2a site Ti'^^ cations. We also found a marginal in- 
crease in the lattice constants with a = b = 5.7216A and 
c = 14.0768A. Susceptibility for the QN sample mea- 
sured with field 5 kOe in the temperature range 2 — 400 K 
(Fig. [2]) shows no sign of magnetic ordering. Data fit- 
ted to CW law in the temperature range (150 — 400 K) 
yields xo = 1.68 x lO"'' cm^/mol, C = 0.145cm3 K/mol 
and 9cw = — lllK. The 9cw is somewhat smaller than 
in the AP sample while C is nearly unchanged. Inverse 
susceptibilities (after subtracting the temperature inde- 
pendent part Xo) of AP and QN samples are linear in 
temperature and deviate from linearity below ^ 80 K, as 
shown on the right axis of Fig. [2l 

The large, negative 9cw indicates that Ir^+ magnetic 
moments are strongly antiferromagnetically coupled with 
each other. Apparently something prevents long range 
magnetic ordering to set in even at 0.35 K (evident from 
heat capacity measurement) which is nearly four hundred 
times lower than 9cw ■ This suggests that inspite of the 
depletion of magnetic ions from the triangular planes, 
geometrical frustration continues to exist in the depleted 
triangular lattice and plays a dominant role in determin- 
ing the magnetic properties of this system. Note that a 
part of the Curie term could be arising from a few per- 
cent of uncorrelated Ir"*"*" spins (and possibly some Ti'^+ 
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Figure 3: Heat capacity of AP and QN Ba3lrTi209 sam- 
ple measured in various applied magnetic fields are shown 
in a log-log scale. Inset: (a) Solid circles represent [Cp(9T)- 
Cp(0 T)]/T of the AP sample and the solid line is the fit (see 
the text), (b) A/fcs as a function of /ioH from OT to 9T is 
shown and the solid line is a fit to Zeeman splitting. 



as well) present in the system, which we call orphan spins 
(discussed later). 

One should note that in literature | 0cw | and ^J■eff 
reported for BasIrTi^Og are greater than 400 K and 
1.73/zb, respectively [l3| which are at variance from our 
data. To clarify this discrepancy, we have fitted the pub- 
lished data (Fig. [9] in supplemental material [3) with 
CW law and found xob = 3.42 x 10"" cmVmol, C = 
O.lOcm^K/mol (^e// = 0.89^*5) and 9cw = -104 K. 
Apparently, Bryne et al. used xoA = 0.5 x 10~^cm'^/mol 
leading them to infer a different ^cw and fief f (see sup- 
plemental material [3| for details). 

Next, in Fig. [3] we present the heat capacity (Cp) 
in various fields for the AP and QN samples (data for 
all fields are shown in supplemental material [lH). No 
anomaly indicative of long range ordering is found in the 
measurement range (0.35 — 295 K). For both the samples, 
Cp depends on the applied field below ~ 20K. This field 
dependence could be arising from a Schottky anomaly of 
orphan spins. We model the heat capacity of Ba3lrTi209 
as arising out of four contributions. These are namely, 
the magnetic contribution of the correlated spins (Cm), 
the lattice contribution (Cat) and the Schottky anomaly 
of the orphan spins (Cgch-orp) and the nuclear Schottky 
anomaly. 

To extract the magnetic part of the heat capac- 
ity arising from correlated magnetic moments, we pro- 
ceed as follows. The Cp has contributions from Cm, 
Ciat, the Schottky anomaly (Csch-orp) from Ir orphan 
spins and nuclear Schottky anomaly (Csch-nuc)- Cat 
is field independent while the others might be field de- 
pendent. Using the zero field heat capacity [Cp(OT)] 
and that measured with 'nT' field [Cp(nT)], we obtain 
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Figure 4: Magnetic heat capacity for the AP and the QN sam- 
ples are shown. The solid lines are fit to power law with power 
indicated in the figure. In the low temperature region, solid 
lines with similar color are with same power. Inset: Magnetic 
entropy change ASm is shown as a function of temperature 
for OT. 



ACp_i,./r = [Cp(iiT)-Cp(OT)]/T. This is then fitted 
with /[Csch(Ai) - Csch(A2)]/T, where / is the percent- 
age of orphan spins present in the sample. Csch(Ai), 
and Csch(A2) are the Schottky anomalies from S = 1/2 
spins and Ai and A2 are the level splittings with applied 
magnetic fields Hi and iJ2, respectively. Here, 
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where R is the universal gas constant and ks is the Boltz- 
man constant. Inset (a) of Fig. [3] shows ACp_ir/T ob- 
tained for OT and 9T along with the fit described above. 
The good fit above ^ 2 K suggests that Cm is not field 
dependent at least above ~ 2 K and all the field depen- 
dence is in Csch-orp- However, below ~ 2K, there is de- 
viation of the fit from the data (this is much larger than 
the expected nuclear Schottky anomaly) which suggests 
the Cm might be field dependent there. The fraction 
of orphan spins / is found to be ~ 3%. The Schot- 
tky splitting (A/ks) obtained from fitting for differ- 
ent fields is plotted as a function of field in the inset 
(b) of Fig. [3l Similar analysis has been reported for 
Ba3CuSb209 i, ZnCu3(OH)6Cl2 M, Y2BaNi05 HH 
etc. At zero field also we found a level splitting of 1.8 K 
which is unexpected but found in Ba3CuSb2 09 @ as 
well. For ^qH > 2T, the Schottky splitting gap follows 
A = gfj^BH, as expected for free spin Schottky anomalies. 
The 'g' value for orphan spins obtained from the linear 
fit is 2.06. 

Using Eq. [TJ the Schottky heat capacity can now 
be subtracted from the measured heat capacity of 



Ba3lrTi20g. Next, we would like to extract the lattice 
heat capacity and for that we have used Ba3ZnSb20g 
as non-magnetic analog. Since the Debye frequency 
is primarily determined by the lighter atoms (in these 
cases oxygens), it will not vary much between these two. 
The high-temperature heat capacities of Ba3lrTi209 and 
Ba3ZnSb2 09 differ because of the difference in their 
molecular weights and lattice volume. By scaling the 
heat capacity of Ba3ZnSb209 (obtained from Ref. [^) 
by a factor of ~ 0.75 we find that the heat capacities 
of Ba3ZnSb209 and Ba3lrTi209 match in the temper- 
ature region ^ 20K-30K. The scaled heat capacity of 
Ba3ZnSb209 is then subtracted from that of Ba3lrTi209 
in order to obtain the magnetic heat capacity as shown 
in Fig. m 

The Cm for both AP and QN samples is independent 
of field from ^ 2.5-10 K and in this range they follow a 
power law in temperature with power 0.4 for the AP 
sample and ^ 0.7 for the QN sample. Above ~ lOK 
the results can be largely affected by uncertainties asso- 
ciated with the subtraction process. Notably, Cm for the 
QN sample is larger than that for the AP sample. Below 
~ 2K, Cm becomes field dependent (for both samples) 
but follows a power law with temperature with the same 
power for different fields. This power is 1.9 for the AP 
sample and 6.5 for the QN sample at very low temper- 
ature (shown in Fig. [5]). From the heat capacity data 
of Ba3CuSb209 (with space group P63/TOmc) published 
in Ref. |22j, we have extracted Cm by subtracting a 
Schottky contribution. Here also we found Cm to be 
field dependent below ^ 2K but following a power law 
with power 2.1 for different fields and field independent 
in the rage 5 — 15K (see supplemental material [16]). In 
many other frustrated systems Cm follows a power law 
with temperature. The power is 2 for the 2D 5 = 1 
system NiGa2S4 0, 1 and 2 for Ba3NiSb2 09 6H-B and 
3C phases respectively, between 2 and 3 for Na4lr308 
Q and 1 at low-temperature but 2 at higher tempera- 
ture in S" = 1/2 system Ba3CuSb209 (with space group 
Pfemc) 01. A power of 2/3 was predicted by Motrunich 
|23| for S — 1/2 triangular lattice organic spin liquid 
system K-(ET)2Cu2(CN)3. In view of the fact that our 
Ir-based system is expected to have a significant spin- 
orbit coupling, a fresh theoretical effort in this direction 
is warranted. 

Magnetic entropy change (A5m) is obtained by inte- 
gration of Cm/T with T and is shown as a function of 
temperature in the inset of Fig. S) The ASm is an or- 
der of magnitude lower than i?ln2 expected for ordered 
S ~ 1/2 systems. In many geometrically frustrated sys- 
tems it is observed that the entropy change is lower 
than expected value. For example, the entropy change 
is 30% and 41% of Mn(2S' + 1) for Ba3CuSb209 and 
Ba3NiSb209 (6H-B phase) respectively, which are similar 
in structure with Ba3lrTi209. However, in Ba3lrTi209 
the magnetic moments are strongly reduced probably due 
to a strong spin-orbit coupling. Here S is not a good 
quantum number and probably Jeff is. So the expected 
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entropy change may not be i?ln(2S' + 1) i.e i?ln2, but 
rather a much smaller quantity. Interestingly, the heat 
capacity is different for the QN sample compared to the 
AP sample implying the influence of atomic site disor- 
der on the details of the triangular lattice and hence the 
ground state. 

Conclusions: We have presented a potentially new 
spin liquid system Ba3lrTi2 09 which is based on a tri- 
angular lattice of Ir^+ ions with electrons responsible for 
the magnetic properties coming from the 5d electronic 
orbital. The sample crystallizes in PGamc space group 
with a large disorder between Ti^+ and Ir^+ cations re- 
sulting in a site dilution of nearly 1/3 of the Ir^"'' sites of 
the edge-shared triangular plane by non-magnetic Ti^+. 
Apparently, magnetic correlations and frustrations are 
still maintained with the absence of magnetic ordering 
down to 0.35 K inspite of a high 6cw value — 130K). 
Associated with this is a magnetic heat capacity which, 
though field dependent, follows a power law with power 
1.9 in the low-temperature range. The QN sample has a 



different behavior. This is somewhat like in Ba3CuSb20g 
where different atomic arrangements (Ref. Q and Ref. 
[2^ ) give rise to different magnetic heat capacity. As 
Nakatsuji et al. [2^ has reported that due to site shar- 
ing between Cu^"'" and Sb^+ ions, a distorted honeycomb 
lattice is formed in Ba3CuSb209, we speculate that a 
similar situation may occur in Ba3lrTi209 yet maintain- 
ing a spin liquid ground state. With the demonstration of 
the existence of a Jg^ y = 1/2 state (having a large spin- 
orbit coupling) in Sr2lr04 [l^l, Ba3lrTi209 is possibly 
an example of a Je// = 1/2 quantum spin liquid system 
and a analog of Ba3CuSb209. This should open up 
a new area pertinent to the search for exotic magnetic 
behaviour in 5d transition metal based compounds. 

Acknowledgement: We thank Department of Sci- 
ence and Technology, Govt, of India for financial sup- 
port. FCC acknowledges the support from National Sci- 
ence Council of Taiwan under project number NSC-100- 
2119-M-002-021. 



[1] P.W. Anderson, Mater. Res. Bull. 8, 153 (1973) 
[2] L. Balents, Nature (London) 464, 199 (2010) 
[3] For a review, see Introduction to Frustrated Magnetism, 
edited by C. Lacroix, P. Mendels, F. Mila, Springer, Hei- 
delberg (2010) 

[4] A. P. Ramirez, B. Hessen, and M.Winklemann, Phys. 

Rev. Lett. 84, 2957 (2000) 
[5] J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. 

M. Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, 

J.-H. Chung, D. G. Nocera, and Y. S. Lee, Phys. Rev. 

Lett. 98, 107204 (2007) 
[6] S. Nakatsuji, Y. Nambu, H. Tonomura, O. Sakai, S. 

Jonas, C. Broholm, H. Tsunetsugu, Y. Qiu, and Y. 

Maeno, Science 309, 1697 (2005) 
[7] Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and 

G. Saito, Phys. Rev. Lett. 91, 107001 (2003) 
[8] Y. Okamoto, M. Nohara, H. Aruga-Katori, and H. Tak- 

agi, Phys. Rev. Lett. 99, 137207 (2007) 
[9] H. D. Zhou, E. S. Choi, G. Li, L. Balicas, C. R. Wiebe, 

Y. Qiu, J. R. D. Copley, and ,1. S. Gardner, Phys. Rev. 

Lett. 106, 147204 (2011) 
[10] J. G. Cheng, G. Li, L. Balicas, ,1. S. Zhou, J. B. Goode- 

nough, Cenke Xu, and H. D. Zhou, Phys. Rev. Lett. 107, 

197204 (2011) 

[11] B. J. Kim, H. Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, C. 

S. Leem, J. Yu, T. W. Noh, C. Kim, S.-J. Oh, J.-H. Park, 

V. Durairaj, G. Cao, and E. Rotenberg, Phys. Rev. Lett. 

101, 076402 (2008) 
[12[ H. Okabe, M. Isobe, E. Takayama-Muromachi, A. Koda, 

S. Takeshita, M. Hiraishi, M. Miyazaki, R. Kadono, Y. 

Miyake, and J. Akimitsu, Phys. Rev. B 83, 155118 (2011) 
[13[ R. C. Byrne, and C. W. Moeller, J. Solid State Chem. 2, 

228 (1970) 

[14[ C. Maunders, J. Etheridge, N. Wright, and H. J. Whit- 
field, Acta Crystallogr., Sect. B: Struct. Sci. 61, 154 
(2005) 

[15[ T. Sakamoto, Y. Doi, and Y. Hinatsu, J. Solid State 



Chem. 179, 2595 (2006) 
[16[ See Supplemental Material below 

[17| G. Radtke, C. Maunders, A. Saul, S. Lazar, H. J. Whit- 
field, J. Etheridge, and G. A. Botton, Phys. Rev. B 81, 
085112 (2010) 

[18[ Y. Singh and P. Gegenwart, Phys. Rev. B 82, 064412 
(2010) 

[19[ S. Chikara, O. Korneta, W. P. Crummett, L. E. DeLong, 
P. Schlottmann, and G. Cao, Phys. Rev. B 80, 140407(R) 
(2009) 

[20[ M. A. de Vries, K. V. Kamenev, W. A. Kockelmann, J. 
Sanchez-Benitez, and A. Harrison, Phys. Rev. Lett. 100, 
157205 (2008) 

[21[ A. P. Ramirez, S.W. Cheong, and M. L. Kaplan, Phys. 
Rev. Lett. 72, 3108 (1994) 

[22[ S. Nakatsuji, K. Kuga, K. Kimura, R. Satake, N. 
Katayama, E. Nishibori, H. Sawa, R. Ishii, M. Hagiwara, 
F. Bridges, T. U. Ito, W. Higemoto, Y. Karaki, M. Halim, 
A. A. Nugroho, J. A. Rodriguez-Rivera, M. A. Green, C. 
Broholm, Science 336, 559 (2012) 

[23[ O. I. Motrunich, Phys. Rev. B 72, 045105 (2005) 

[24| B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, 
H. Takagi, and T. Arima, Science 323, 1329 (2009) 



5 



Supplemental material for "Spin liquid 
behaviour in Jgjj = 1/2 triangular lattice 
BaglrTisOg" 

The XRD refinement of the AP sample is sliown in Fig. 
[5l The crystal parameters obtained from refinement are 
given in Table H] We have studied the change in the re- 
finement parameters by varying the site disorder of Ir'*^ 
ions with Ti^'^ ions at Ti(2) and Ti(3) sites. The param- 
eters thus obtained are shown in Fig. [6]and Fig. [71 From 
the figures we can conclude that (37± 10)% disorder with 
Ti(2) site and (7±4)% disorder with Ti(3) site gives best 
refinement. 

ZFC and FC susceptibilities of the AP sample are 
shown for different fields in Fig. |S1 

As mentioned in the main paper, the fJ-eff and | ^cw I 
values obtained by Bryne et al. [l3| are much higher 
than found in our measurements. To find the reason of 
this mismatch, we have reanalysed the published data on 
Ba3lrTi20g in Ref. [l^. The blue open squares in Fig. 
ini represent the inverse susceptibility data as published 
in Ref. (l3| . The corresponding susceptibility is shown 
on the left axis as blue solid squares. We have fitted 
this susceptibility data with the Curie- Weiss (CW) law 
in the range 77 — 363 K. This fitting yields temperature 
independent susceptibility xqb = 3.42 x 10"'' cm'^/niol, 
Curie constant C = 0.10 cm^^K/mol (^fe// — O.SO/is) and 
Curie- Weiss temperature 6cw = — 104 K which are some- 
what closer to the values obtained from our measure- 
ment. Subtracting this xob, we have plotted (% — Xob)~^ 
as pink diamonds which is linear in the whole tem- 
perature range. Apparently, Bryne et al. did not fit 
their susceptibility data with CW law to find C and xo, 
rather they have chosen a temperature independent sus- 
ceptibility xoA = 0.50 X 10~*cm^/mol. We have also 
shown (x — Xoa)~^ data points as green solid circles 
which are much different from (x — Xob)^^ data points. 
The slope corresponding to these green circles has been 
used by Bryne et al. jl3| to find ^cw- Further they 
have used this 6cw and a single (at 293K) susceptibil- 
ity data point (xm) S^^ Me// from the formula 
fieff = 2.83-y/xM(r - 9). Hence, fixing xoA (without a 
fitting procedure) and calculating /ig/ / based on a single 
susceptibility data point gives unreliable /ie/ / and ^cw 
values in Ref. [T3| . 

Heat capacities for different fields for the AP and QN 
sample are shown Fig. [10] and Fig. [11] respectively. 
Schottky fits for different fields for the AP sample are 
shown in Fig. \T% 

In a recent report on Ba3CuSb20Q (space group 
PGs/mmc), Nakatsuji et al. (Ref. [22|) have shown 
the heat capacity of the system after subtracting the 
lattice contribution. We term this as C^. We have 
analysed their data to extract magnetic heat capacity 
(Cm) after subtracting the Schottky contribution. Inset 
of Fig. [131 shows [C^^(5T)-C^(0T)]/r and its fit with 
/[Csch(Ai) — Csch(A2)]/r (see main paper for details). 
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Figure 5: Refinement pattern of as-prepared (AP) 
Ba3lrTi209 is shown. Inset: Normalised main peaks for the 
AP and the quenched sample (QN) are shown. 



Table I: Atomic parameters obtained by refining x-ray powder 
diffraction for as-prepared Ba3lrTi2 09 at room temperature 
with a space group PGsmc. 





X 


y 


z 


g 


Ba(l) 


2a 








0.24199 


1.00 


Ba(2) 


2b 


1/3 


2/3 


0.07981 


1.00 


Ba(3) 


2b 


1/3 


2/3 


0.39153 


1.00 


Ir(l) 


2b 


1/3 


2/3 


0.64464 


0.56 


Ti(l) 


2b 


1/3 


2/3 


0.64464 


0.42 


Ti(2) 


2b 


1/3 


2/3 


0.83775 


0.63 


Ir(2) 


2b 


1/3 


2/3 


0.83775 


0.37 


Ti(3) 


2a 








0.48763 


0.93 


Ir(3) 


2a 








0.48763 


0.07 


0(1) 


6c 


0.16098 


0.83898 


0.57538 


1.00 


0(2) 


6c 


0.48859 


0.51138 


0.74989 


1.00 


0(3) 


6c 


0.16098 


0.83898 


0.91548 


1.00 



The fit is good above ^ 1.5 K which means Cm is field 
independent and the field dependence of C^ in that re- 
gion is totally coming from Schottky contribution. Below 
^ 1.5 K, the fit deviates from data points indicating field 
dependence of Cm ■ The A obtained from the fit has been 
used in Eq. 1 of our paper to get the Schottky contribu- 
tion to heat capacity (/ was inferred to be 17%). This 
is subtracted from C^ to get Cm for Ba3CuSb209. For 
different fields we have extracted Cm as shown in Fig. 
[T31 Cm is field independent in the range 5 — 15K and 
follows a power law with temperature with power 1. Be- 
low 2 K Cm is field dependent but still follows a power 
law with power 2.1 for different fields. The behavior is 
very similar to that found by us in Ba3lrTi209. 
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Figure 6: Refinement parameters for the as-prepared (AP) 
sample obtained by varying disorder at Ti(2) site keeping dis- 
order at Ti(3 )site unciianged is siiown. Rp and (Rmp-2) corre- 
sponds to the left axis while goodness of fit (GOF) is plotted 
on right been offset downward by 2 to show Rp 

and Kwp on the same axis with clarity. 
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Figure 7: Refinement parameters for the as-prepared (AP) 
sample obtained by varying disorder at Ti(3) site keeping dis- 
order at Ti(2) site unchanged is shown. Rp and (Rtup-2) corre- 
sponds to the left axis while goodness of fit (GOF) is plotted 
on right axis. R^jp has been offset downward by 2 to show Rp 
and Kwp on the same axis with clarity. 
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Figure 8: Magnetic susceptibilities of as prepared (AP) 
Ba3lrTi209 sample measured at different fields are shown in 
a semilog scale. 




Figure 9: Inverse susceptibility data of Ba3lrTi2 09 obtained 
from Ref. is shown as blue open squares. The corre- 

sponding susceptibility is shown as blue solid squares with its 
fit with Curie Weiss (CW) law (green solid line). Inverse sus- 
ceptibility after subtracting xoa = 0.5 x 10"'* cm^/mol (as in 
Ref. [31) is shown as green solid circles. The inverse suscepti- 
bility after subtracting xob = 3.42 x 10~*cm^/mol (obtained 
from CW fit done by us) is also shown as pink diamonds. The 
dashed line is a guide to eye. Susceptibility is plotted on the 
left axis while the inverse susceptibilities correspond to the 
right axis. 
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Figure 10: Heat capacities of as prepared (AP) Ba3lrTi209 
sample measured in various applied magnetic fields are shown 
in a log-log scale. 




Figure 11: Heat capacities of quenched (QN) Ba3lrTi209 
sample measured in various applied magnetic fields are shown 
in a log-log scale. 
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Figure 12: Scattered symbols represent [Cp(nT)-Cp(OT)]/r 
of the AP sample and the solid line of corresponding color is 
the fit (described in the main paper). 




Figure 13: Magnetic heat capacities (Cm) of Ba3CuSb209 
(taken from Ref. (2^ ) for different fields after subtracting 
Schottky contribution are shown. The solid lines are fits to 
power law. In the low-temperature region Cm is field depen- 
dent and the power is 2.1 for different fields. In the range 
5 — 15 K, Cm is independent of field and linear in temper- 
ature. Inset: Solid squares represent [Cm(5 T)-Cm(0 T)]/r 
and the solid line is the fit (see text). 



10 



